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Abstract 

In total 1317 di-lepton events with significant missing transverse momentum were iden- 
tified in a total data sample of 680 pb~^ collected at e"'"e~ centre-of-mass energies ranging 
from 183 GeV to 209 GeV. The number of di-lepton events, the dependence on centre- 
of-mass energy, and the event properties are consistent with expectations from Standard 
Model processes, predominantly W"'"W~ production with both W bosons decaying lep- 
tonically. This topology is also an experimental signature for the pair production of new 
particles that decay to a charged lepton accompanied by one or more invisible particles. 
No evidence for new phenomena is apparent. Upper limits are presented on the production 
cross-section multiplied by the relevant branching ratio squared for sleptons, leptonically 
decaying charginos and charged Higgs bosons. Mass limits are also given. 
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1 Introduction 

We report the final results from an investigation of events containing two oppositely charged 
leptons and significant missing transverse momentum recorded with the OPAL detector at 
LEP. These events were produced between 1997 and 2000 in the highest energy e+e^ collisions 
achieved to date, at centre-of-mass energies {^/s) ranging from 183 to 209 GeV. The data 
used in this analysis amount to a total integrated luminosity of 680 pb~^. The event topology 
which is studied, henceforth called acoplanar di-lepton events, consists of low multiplicity events 
with two oppositely charged leptons, significant missing transverse momentum and the possible 
presence of additional photons. 

The number of observed events and their properties are compared with the expectations 
for Standard Model processes, which are dominated by the i'^i/i'V final state {i = e, /x, r) 
arising from W+W~ production in which both W bosons decay leptonically. This topology is 
also an experimental signature for the pair production of new particles that decay to produce 
a charged lepton accompanied by one or more invisible particles. The invisible particles may 
be neutrinos, or the lightest stable supersymmetric [1] particle (LSP), or weakly interacting 
neutral particles with long lifetimes, which decay outside the detector volume. The LSP may 
be the lightest neutralino, Xi: the lightest sneutrino, i>, or the gravitino, G. We present the 
results of searches for the pair production and stated decay mode of the following new particles: 

charged scalar leptons (sleptons): i"^ ^^Xi (or ~* ^^G), where may be a selectron 
(e), smuon (/i) or stau (f) and £^ is the corresponding charged lepton. 

chcirged Higgs: H='= t^Ut-. 

chcirginos: xf ~^ ("2-body" decays) or xf — > i^vxi ("3-body" decays). 

The slepton searches are also relevant to interpreting the results of searches for chargino 
and neutralino production since the chargino and neutralino production cross-sections and 
branching ratios depend on the sneutrino and charged slepton masses. The search for charged 
sleptons provides constraints on the slepton masses, notably the selectron mass, while indirect 
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limits on the sneutrino masses, notably the electron-sneutrino mass, can be obtained in models 
where the charged slepton and sneutrino masses are related. The search results described here 
will also be used in a separate publication regarding the search for sleptons with non-negligible 
lifetime. 

In most respects the analysis is similar to our published searches at centre-of-mass energies 
of 161, 172, 183 and 189 GeV [2-4]. This paper supersedes the results of [3,4]. The analysis 
is performed in two stages. The first stage consists of a general event selection designed for 
all possible low multiplicity events containing a lepton pair plus missing transverse momentum 
(Section 3). In this context the Standard Model C.^vC~V events are considered as signal in 
addition to possible new physics sources. All Standard Model processes that do not lead 
to i^ui'u final states, e.g. e~'"e~£+£" and £"{•-)') , are considered as background and are 
reduced to a rather low level (r^ 3%) by the event selection. The observed numbers of events 
and kinematic distributions for the data are compared with expectations from Standard Model 
processes in Section 4. In the second stage the detailed properties of the events (e.g. the 
type of leptons observed and their momenta), which vary greatly depending on the type of 
new particles considered and on free parameters within the models, are used to separate as far 
as possible the events consistent with potential new physics sources from W''"W~ production 
and other Standard Model processes (Section 5). Constraints on new physics are discussed in 
Section 6. 

Slepton search results from the ALEPH collaboration at ^/s <209 GeV [5] have been pub- 
lished recently. The L3 and DELPHI collaborations have published searches for sleptons in this 
channel using data with ^/s <189 GeV [6]. 

2 OPAL Detector and Monte Carlo Simulation 

A detailed description of the OPAL detector can be found elsewhere [7, 8] and only the general 
features are described here. 

The central detector consisted of a system of chambers providing charged particle tracking 
over 96% of the full solid angle inside a 0.435 T uniform magnetic field parallel to the beam 
axis. It consisted of a two-layer silicon micro-strip vertex detector, a high precision vertex 
drift chamber, a large volume jet chamber and a set of 2;-chambers that measured the track 
coordinates along the beam direction. 

A lead-glass electromagnetic calorimeter (ECAL) was located outside the magnet coil. It 
provided, in combination with the gamma-catchers (GC) and forward detectors (FD), which 
were lead-scintillator sandwich calorimeters and, at smaller angles, silicon tungsten calorimeters 
(SW), geometrical acceptance with excellent hermeticity down to approximately 25 mrad^ from 
the beam direction. 

The magnet return yoke was instrumented for hadron calorimetry and consisted of barrel 
and endcap sections along with pole tip detectors that together covered the region | cos 6'| < 0.99. 
Outside the hadron calorimeter (HCAL), four layers of muon chambers covered the polar angle 
range of |eos^| < 0.98. Arrays of thin scintillating tiles were installed in the endcap region 

^For some polar angles, precision energy measurements were compromised by upstream material such as the 
synchrotron shield installed for LEP2 at approximately 30 mrad. However, the ability to veto significant energy 
deposits with very high efficiency extended down to polar angles of 25 mrad. 
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to improve trigger performance, time resolution and hermeticity for operation at LEP 2 [8]. 
Of particular relevance to this analysis are the four layers of scintillating tiles (the MlP-plug) 
installed at each end of the OPAL detector covering the angular range 43 < ^ < 220 mrad. 
These tiles were commissioned in 1997 and became fully operational for data taken from 1998 
at ^/s > 189 GeV. Time-of-flight (TOF) scintillators in the barrel region aided cosmic ray 
rejection. 

The following Standard Model processes are simulated. Four-fermion production is simu- 
lated using the KORALW [9] generator which uses the grc4f [10] matrix elements to calculate 
the four-fermion cross-sections including interference effects and includes a detailed descrip- 
tion of hard radiation from initial and final state charged particles. Two-photon production of 
muon pairs and tau pairs is generated with the BDK generator [11] and the program of Ver- 
maseren [12] is used for the two-photon production of electron pairs. The grc4f event generator 
excluding multi-peripheral processes is used for e+c^f+f^and e+e^qq. The Phojet [13] and 
Herwig [14] event generators are used to study backgrounds from two-photon production of 
hadrons. The production of lepton pairs is generated using Bhwide [15] and Teegg [16] for 
e"'"e~(7), and KK2F [17] for and t~^t~{j). The final state 1/1/7(7) is generated with 

NUNUGPV98 [18]. 

Slepton pair production is generated using Susygen [19]. Chargino pair production is 
generated using Dfgt [20] for three-body decays, and Susygen for two-body decays. Charged 
Higgs boson pair production is generated using HzHA [21]. 

Backgrounds from the accelerator or cosmic-ray interactions can lead to additional hits, 
energy deposition and even reconstructed tracks being superimposed on triggered data events. 
Such effects of detector occupancy are simulated for all Monte Carlo samples by adding to 
the Monte Carlo events the hits, energy depositions and additional jets found in randomly 
triggered beam-crossing data events (BXRSA trigger [22] with a constant 0.1 Hz rate) corre- 
sponding to the same centre-of-mass energy. Systematic checks showed that this constant rate 
model adequately described the average inefficiency during a LEP fill, despite the instantaneous 
luminosity decreasing typically by a factor of 2 or 3. 

All Standard Model and new physics Monte Carlo samples are processed with a full simula- 
tion of the OPAL detector [23] and subjected to the same reconstruction and analysis programs 
as used for the OPAL data. 



3 General Event Selection 

3.1 Introduction 

The general event selection is formed by requiring that an event be selected by either or both 
of two independent event selections, referred to here as selection I and selection IL Selection 
I is designed to retain efficiency for events with low visible energy. This is characteristic of 
slepton or chargino events where the mass difference. Am , between the parent particle and 
the invisible daughter particle (e.g. Xi) is small. Selection 11 is optimised to maximise the 
efficiency for W+W"— > £+i/£~I7 events, while keeping other Standard Model background events 
to a minimum. 
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Both selections require evidence for two^ charged leptons and an invisible system carrying 
significant missing transverse momentum (p™'^^). The maximum p"^^^^ which can be carried 
away by undetected particles travelling close to the beam is set by the maximum angle to the 
beam at which such a particle will not be detected. This is 25 mrad - the angle to which 
the silicon tungsten detector extends. A particle with energy of E'beam may thus carry away 
puuss _ 0.025 -Bbeam without being detected. 

Some background processes containing secondary neutrinos (particularly from tau decay) 
may have large values of p™'^^ with the direction of the missing momentum vector pointing 
away from the beam axis. Such events tend to be fairly coplanar, and the component of p™'**** 
which is perpendicular to the event thrust axis in the transverse plane (called a"^^^^) is much 
less sensitive than p™'^^ to the presence of neutrinos from tau decays or to poorly measured 
particles. This can be seen by considering electrons produced in tau decay. Low momentum 
electrons produced from this source can have a large angle relative to the original tau direction, 
but their momentum transverse to the original direction (and hence their contribution to a^^^) 
is small. 

The event selection as it was first implemented was described in detail in [2]. In [3] we made 
use of the improved hermeticity for non-showering particles in the forward direction provided 
by the MlP-plug^. Subsequent improvements have been made for the analysis of the data 
taken at 189 GeV [4] by reducing the sensitivity to mis- measurements of p™^^^ and using the 
much improved background rejection capabilities of the MlP-plug to increase substantially the 
efficiency of selection 11. For this paper, selection 1 has been consolidated and much refined. In 
particular, a selection has been added for single lepton events in which a single high transverse 
momentum charged lepton is observed in the central detector. These are usually events where 
the second charged lepton is produced at a polar angle sufficiently close to the beam direction 
that it does not produce a visible track in the central detector. In the context of this paper, this 
"single-lepton" selection is used to recuperate di-leptons where the second lepton was missed. 
Various small modifications were made to selection II in order to reduce sensitivity to poorly 
modelled backgrounds. An example of this is using stricter track quality criteria to reduce 
the possibility of tracks from overlaid beam-gas events affecting the measurement of the event 
kinematics. 

The current event selection is summarised below. 
3.2 Selection I 

Selection I is designed to retain efficiency for events with very low visible energy, but nevertheless 
significant p^^^^^. This is typical of new physics signal events with small Am. The selection 
requires evidence that a pair of leptons has been produced and of significant missing transverse 
momentum. Subsequent cuts reduce the probability that the signature of missing transverse 
momentum is faked by events with secondary neutrinos from tau decay or poorly measured 
particles. 

^This is not strictly true for the single lepton selection of selection I (Section 3.2.3). 

^Note that the general event selection and the Monte Carlo simulation samples used for the 183 GeV data- 
set axe unchanged from the 183 GeV publication. The poorer MlP-plug detector performance in 1997 would 
have necessitated a re-optimisation of the latest general event selection, and we judged such a re-analysis to be 
inappropriate. 



6 



At least one lepton in the event is required to be well identified and to satisfy requirements 
on isolation and transverse momentum. Much looser requirements are made on the possible 
presence of a second lepton in the event. 

In order to be considered in the event selection, tracks in the central detector are required 
to satisfy pt > 0.1 GcV. Clusters in the barrel ECAL are required to satisfy -E > 0.1 GeV 
and clusters in the endcap ECAL are required to satisfy E > 0.2 GeV. Converting photons 
are identified and the tracks and clusters associated to the conversion are replaced by a single 
4- vector representing the photon. 

3.2.1 Lepton Candidates 

The first stage is to look for lepton candidates in the event. A track is identified as a lepton 
candidate if it has p > 1.5 GeV and it is identified as an electron, muon or hadronic tau decay"^. 
The electron ID is based on the ratio of ECAL energy to track momentum (E/p), and dE/dx 
information. Muons are identified using muon chamber or HCAL hits which match with a 
track in the central detector, or from a high momentum track which matches with a low energy 
ECAL cluster. To identify a hadronic tau the following criteria are applied: 

• Within a cone of half-opening angle 35° around a track oi pt >1.5 GeV there are no more 
than three tracks in total. 

• The invariant mass of all tracks and clusters within the cone is less than the tau mass. This 
is calculated assuming the pion mass for each track and correcting for double counting of 
tracks and clusters. 

The lepton candidates are also required to be isolated. There must be no more than 2 
additional tracks and no more than 2 additional clusters in an isolation cone defined around 
the lepton candidate (half opening angle 20° for electrons and muons, half opening angle 60° 
for hadronic tau decays) , and the energy sum of the tracks or of the clusters must be less than 
2 GeV. 

3.2.2 Di-Iepton Event Selection 

Firstly evidence for a pair of leptons is required: 

• There must be at least one and no more than two isolated lepton candidates with pt > 
1.5 GeV. 

• If the event contains two isolated lepton candidates then all tracks in the event must be 
associated with the lepton candidates. 

• If there is only one isolated lepton candidate then the other tracks and clusters in the 
event are considered as a possible second lepton candidate provided the following: 

"^The lepton identification (ID) made at the event selection stage is used for the purpose of the selection only. 
A separate lepton ID is applied to selected events and used by the search analysis (see section 4.1). 
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(a) There must be between 1 and 3 additional tracks, at least one of which must have 
Pt > 0.3 GeV. 

(b) The invariant mass of the additional tracks must be less than 3 GeV and the invariant 
mass of the additional tracks and clusters must be less than 8 GcV. 

(c) 7/9, the net momentum of the additional tracks and clusters divided by their invariant 
mass, is required to be greater than 2.0. 

Next, significant missing transverse momentum is required. For events with a large acopla- 
narity angle (0acop > 7r/2 where 0acop is the supplement of the azimuthal opening angle), the 
following cut is applied ^: 

• Xt > 0.045, where Xt is the scaled missing transverse momentum of the event {p^^^^ / E^eam) ■ 
For events with small acoplanarity angle (0acop < ^f^)'- 

• xt> 0.035 

• At small acoplanarity, cuts on a^^'Y-^beam and 6'™*''' (where 6'™*'*' = tan~^[a™''YP™'''1' 
p™'®® is the total momentum of the observed particles in the z direction) are used to reduce 
background from processes such as t'^t~ and e'^e" t'^t~ . Events are divided into subsets 
using variables which help estimate how likely they are to originate from r^r~ or c^c~ 
T~^T~ . The cut values vary according to the subset (in increasing order of similarity to 
r+T" or e+e" r+r" background): 

(a) Xt > 0.2 and unassociated energy less than 1.5 GeV and the lepton pair is identified 
as e+e- , fx+ii' or e^^u^ : af^'/Ei,^^^ > 0.011 and 9f'' > 0.025 

(b) Xt > 0.2 or unassociated energy less than 1.5 GeV but not in subset (a) : a™****/£^beam > 
0.018 and Of'' > 0.05 

(c) Xt < 0.2 and unassociated energy greater than 1.5 GeV: a™^'^'^/-Ebeam > 0.025 and 
9f^^ > 0.1 

• Xt + arV^beam > 0.07. 

Further cuts are apphed to reduce the effect of processes which may fake the signature of 
missing transverse momentum. These are mostly vetoes against energy in the forward region 
(GC, FD, SW or MlP-plug). However care is taken to ensure that the activity in the forward 
region is only used to veto the event if the activity could possibly explain the apparent miss- 
ing momentum. Also different requirements are made depending on the amount of missing 
transverse momentum observed. 

Radiative lepton pair events containing high energy isolated photons form a potentially 
serious source of background, because quantities such as p™'^^ and a™'^^ may be poorly measured. 
However, in order to retain efficiency for potential signal events containing isolated photons, 
events are rejected only if the presence of the isolated photon could possibly have caused the 

observed p"^^^^ and af^^^^ . 

^The effect of measurement errors is taken into account in selection I by taking each lepton in turn and 
shifting its momentum up and down by one standard deviation of its estimated measurement error. At each 
stage the values of p™^^'^ and a™'"" are recalculated and the minimum of all the values obtained is the one used 
for comparison to the cut values. 
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3.2.3 Single Lepton Selection 

Selection I also includes a selection for events with only one lepton visible in the central detector, 
and the other lepton travelling sufficiently close to the beam axis that no track is produced in 
the central tracking chambers. 

The principal single lepton selection criteria are: 

• The event must contain one and only one identified, isolated lepton candidate, and no 
other tracks. 

• xt> 0.16 

In order to veto events which may fake the missing transverse momentum signature, events 
are rejected if they contain activity in ECAL, GC, FD, SW, muon endcap or MlP-plug which 
is back-to-back with the observed lepton. 

3.3 Selection II 

Selection II is optimised to select high visible energy events typical of the W+W~— > t^v l~V 
process. A low multiphcity preselection is apphed such that the events contain at least one 
track but no more than 8. Also the sum of the number of tracks plus the number of ECAL 
clusters is required to be less than 16. 

A cone-based jet finding algorithm [24] is applied requiring a minimum jet energy of 2.5 GeV 
and a cone half angle of 20°. Events are required to contain 1, 2 or 3 jets, and a separate 
selection is used for each value of njet- The majority (about 90%) of W+W"— >• i^u events 
have njet = 2. One-jet events are usually those where the decay products from one of the W 
bosons are poorly reconstructed (for example if the lepton is travelling close to the beam pipe) . 
Three-jet events can occur if there is an energetic photon in the event. 

Electron and muon identification, similar to that in selection I, is applied to the most 
energetic track in each jet. Jets not identified as electrons or muons are classified as hadronic 

tau decays. 

The most important cuts for each njet class are summarised below. 

3.3.1 Di-jet Selection 

1- ^acoi > 5°, where 6'acoi is the acolinearity angle between the two jets (defined as the 
supplement of the three-dimensional opening angle). 

2. Xt > 0.05. It is further required that the significance by which xt exceeds 0.05 be greater 
than 1 standard deviation. 

3. For events with 0acop > 7r/2 it is required that the direction of the missing momentum 
satisfy | cos 6'°"'^'^ | < 0.95. For events with 0acop < 7r/2 it is required that a™'^^/£^beam > 
0.022 and that sin 9f^^ > 0.06. 
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Further cuts are made on the quahty of each jet, and for background rejection (mainly for 
events at low Xt ): 

• Events are rejected if there are any tracks which are not associated with either jet. 

• Events with low Xt {xt < 0.15) are rejected if there is evidence of activity in the MIP- 
plug with azimuthal angle within 60° of the missing transverse momentum direction. 
For events with medium Xt (0.15 < Xt < 0.23), the veto uses only the outer MlP-plug 
scintillators. 

3.3.2 Tri-jet Selection 

For events classified as tri-jet, significant missing momentum is required: 

• The sum of the opening angles among the three two-jet pairings should be less than 359°. 

• Xt of the three-jet system should exceed 0.05 with a significance exceeding 1 standard 
deviation. 

Further requirements depend in part on the thrcc-jet event topology in the transverse plane 
characterised by A0™^^. This is calculated by ordering the jets by increasing azimuth, and 
finding the maximum azimuthal di-jet separation angle when rotating clockwise from one jet 
to the next^. 

The most important additional requirements are: 

• For A(j)"^'^^ > 185°, there should be no pair of vertex drift chamber axial tracks with an 
azimuthal opening angle exceeding 165°. These tracks are reconstructed independently 
from standard jet chamber based tracks and this criterion helps to reduce cases where 
the wrong jet chamber left- right ambiguity is chosen in the standard tracking. 

• For A0™^^ < 185° it is explicitly required that two of the jets have associated tracks 
(charged jet) and the third jet has no associated tracks (neutral jet). For the events with 
A0™^ < 180° an axis in the transverse plane is defined using the most energetic charged 
jet. The event is rejected if the transverse momentum of the neutral jet with respect to 
this axis exceeds 80% of the transverse momentum of the less energetic charged jet. This 
cut is effective against t~^t~j events. 

• Events are rejected if there is evidence for a particle passing through the MlP-plug (similar 
to di-jet MlP-plug veto but without the directional requirement). 

3.3.3 Single-jet Selection 

The single-jet selection applies to events where one lepton with high transverse momentum 
is observed at wide angle | cos6'| < 0.82 with evidence for a partially reconstructed lepton at 
small polar angle, or events where the two leptons fall within the same cone (massive mono-jet) 

^This definition results in A(^™'*^ being in the range from 120° to 360°. 
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for I cos ^ I < 0.90. In contrast to the single lepton category of selection I, it is required here 
that there be some activity in the forward region (endcap ECAL, GC, FD, SW, muon endcap, 
MlP-plug) when there is evidence for only one lepton at wide angle. 

Backgrounds from cosmic rays are reduced by requiring in-time TOF hits (within 20 ns) for 
tracks in the barrel region, and that the most energetic track in the jet be associated with hits 
in the silicon micro-vertex detector. 

Three different minimum requirements on Xt are applied (0.20, 0.30 and 0.25) depending 
on whether the event satisfies the massive mono-jet criteria and whether the forward activity 
is solely due to the MlP-plug or not, respectively. 

4 Comparison of Data with Standard Model Processes 

The numbers of events passing the general selection at each centre-of-mass energy in the data 
are compared with the Standard Model Monte Carlo predictions in Table 1. The total number 
of events predicted by the Standard Model is given, together with a breakdown into the contri- 
butions from individual processes. The data have been binned by centre-of-mass energy in bins 
which reflect the predominant centre-of-mass energies where data was collected. The number 
of observed candidates is consistent with the expectation from Standard Model sources, which 
is dominated by the i'^u £~V final state arising mostly from W''"W~ production. 
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0.31±0.04 
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0.6±0.06 


2.1±0.4 
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2.0±0.3 


0.38±0.07 
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243 


249.5±1.4 


241.4±1.2 
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0.64±0.08 


207.9 


9.03 


19 


18.2±0.2 


17.7±0.2 


0.20±0.04 


O.OSiO.Ol 


0.18±0.03 


0.06±0.01 


All 


680.4 


1317 


1303.6±3.3 


1259.8±2.9 


19.6±1.2 


5.4±0.2 


15.2±0.7 


3.4±0.2 



Table 1: Comparison between data and Monte Carlo of the number of events passing the 
general selection in each centre-of-mass energy bin. The total number of events predicted by 
the Standard Model is given, together with a breakdown into the contributions from individual 
processes. The Monte Carlo statistical errors are shown. Also listed for each bin is the mean 
■s/s (in GeV) and the integrated luminosity, L, in pb~^ . 



4.1 Lepton Identification and Kinematic Distributions 

Information about the type of lepton found in the selected events and measurement of the event 
kinematics was used to compare the data with expectations from Standard Model processes. 
Some of these distributions are then used as hkelihood variables to distinguish between Standard 
Model and new physics sources of acoplanar di- lepton events (section 5). We examine the 
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Lepton ID 


True Identity of Lept 


on 


e or r — >• e 


/X or r — > /i 


h 


e 


96.3 % 


0.1% 


4.7 % 




0.1 % 


98.2 % 


3.5 % 


h 


1.2 % 


0.8 % 


86.9 % 


e/h 


2.2 % 


0.1 % 


2.8 % 




0.2 % 


0.8 % 


1.8 % 


X 


0.0 % 


0.0% 


0.4 % 



Table 2: Lepton identification performance, calculated using KoralW four-lepton Monte Carlo 
events at y/s —189-206 GeV. The efficiency is evaluated for leptons generated with \ cos 9\ < 0.95 
and p/Ey,^^^ > 0.02. 

distributions of the likelihood variables and other related kinematics distributions, as measured 
in the data, and compare to expectations from Standard Model processes. The following 
distributions are used as the likelihood variables: 

• Di-lepton identities; 

• Acolinearity of the event, defined as the supplement of the three-dimensional opening 
angle between the two leptons; 

• Momentum asymmetry of the event, — P2I/ (P1+P2)), where pi and p2 are the momenta 
of each lepton; 

• Scaled momentum of each lepton, (p/£'beam) using the best estimate combining tracking 
and calorimetry information; 

• —q cos 9 of each lepton, where q is the particle charge and 9 is the polar angle with respect 
to the electron beam direction. 

The lepton identification algorithm has been improved using techniques developed for tau 
lepton decay analyses such as described in [25] . The efficiencies for correctly identifying leptons 
as electrons, muons or hadronically decaying tau leptons are shown in Table 2 for a sample of 
leptons representative of the geometrical and kinematical acceptance of the general selection. 
Note that we have six separate mutually exclusive lepton identification decisions: electron 
(e), muon hadronically decaying tau (h), electron-hadron ambiguous (e/h), muon-hadron 
ambiguous (/x/h) and "rest-of-event" (x). The ambiguous categories correspond to cases where 
the lepton candidate has properties which do not permit clear identification between the two 
lepton types, while the rest-of-event category arises mostly from tau leptons which fail to pass 
basic lepton identification criteria. 

We show in Figure 1 the observed di-lepton identities and the above kinematic variables 
are shown in Figure 2. Additional kinematic variables, which are not directly used in the 
likelihood, are also shown in Figure 3. Reasonable agreement is observed between the data and 
the Standard Model Monte Carlo in all these distributions as quantified using tests. 
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5 Likelihood Method for New Physics Search 



5.1 Method 

The method used in the second stage of the analysis, in which we distinguish between Standard 
Model and new physics sources of lepton pair events with missing momentum, is essentially 
described in [4]. Discrimination is provided by information on the lepton identification, the 
acolinearity and momentum asymmetry of the event, the scaled momentum and —qcos9 of the 
observed lepton candidates. These variables are combined using a likelihood technique. The 
discriminating quantity used is the relative likelihood, Lr , defined by: 



where Ls is the likehhood of the event being consistent with the signal hypothesis (a particular 
new physics scenario) and Lb is the likelihood of the event being consistent with the background 
hypothesis (Standard Model sources of lepton pair events with missing transverse momentum). 
The distributions depend on Am and on the parent particle mass, m. 

In our previous publications we justified the inclusion of most of these variables and de- 
scribed how they were used [4]. In particular the —qcosO variable is not used in the selectron 
and chargino searches because the t-channel neutralino (selectron search) and sneutrino ex- 
change (chargino search) contributions make this signal distribution model-dependent. 

The new momentum asymmetry variable describes the correlation between the momenta 
of the two leptons and is thus complementary to the two scaled momentum variables which 
are implemented in the likelihood function as independent variables. It is correlated with the 
acolinearity variable, and so the information was included using a two-dimensional probability 
distribution of momentum asymmetry versus acolinearity. Sensitivity studies for the selectron 
and smuon searches justified this addition. However similar studies for the other channels (stau, 
chargino and charged Higgs), where momentum distributions are broader, showed no net gain in 
sensitivity and so for these channels the acolinearity distribution alone was used. Also for these 
channels, we removed the distinction between leptons being identified as electrons and muons 
and merged them into one category since we found that there was no significant difference in 
search sensitivity. 

For the selectron search, we required that at least one lepton be identified as an electron 
and that no lepton in the event be compatible with being a muon (i.e. the other lepton was not 
identified as a muon and it was not muon-hadron ambiguous). Similarly for the smuon search 
we required at least one identified muon and required that no lepton be compatible with being 
an electron. Some additional requirements were used to reject events from these searches if the 
second lepton was identified as a hadronically decaying tau or a "rest-of-event" candidate and 
if it had properties strongly incompatible with being an electron for the selectron search and 
with being a muon for the smuon search. These additional requirements cut the background 
by 21% and 14% while only reducing the signal efficiency by 2.1% and 0.8% for the selectron 
and smuon searches respectively. 

For each search channel, reference histograms were constructed for each of the likelihood 
variables at each point in m and Am for which signal Monte Carlo had been generated (each 
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point corresponds to a table entry in tables 3-8). A smoothing algorithm [26] was applied to the 

histograms to reduce the effects of statistical fluctuations. The reference histograms were then 
used to construct Lr distributions. An interpolation procedure was developed which allowed 
us to interpolate among neighbouring signal Monte Carlo (m, Am , y/s) points to intermediate 
values of (m, Am and a/s). Further details are given in [27]. 

Lr distributions for signal Monte Carlo, Standard Model Monte Carlo and data are shown 
in Figure 4 for the specific example of the analysis for staus with a mass of 80 GeV and a 
stau-neutralino mass difference of 60 GeV. There is considerable variation in the shapes of 
these distributions with m and Am but agreement between data and Standard Model Monte 
Carlo is generally good. 

5.2 Calculation of Cross-section Upper Limits 

For each search channel and for each set of kinematic parameters (m and Am) we wish to test 
the consistency of the data with the sum of background and an additional signal contribution. 
The signal contribution depends on as, the cross-section multiplied by branching ratio to lep- 
tonic final states squared. We use an extended maximum likelihood technique as described in [4] 
to measure cTg and, in the absence of a significant signal, to set 95% confidence level (CL) upper 
limits on cTg. We form a likelihood, L{as), of the set of Lr values for the data being consistent 
with the expected Lr distribution for Standard Model plus a signal contribution of Ug. The 
upper limit on the cross-section multiplied by branching ratio squared at 95% confidence level, 
(T95, is calculated as the value of as below which 95% of the area under the likelihood function 
lies. Details of the likelihood function and method as well as cross-checks of the method are 
described in [4]. Data from different centre-of-mass energies are combined by weighting them 
using an assumed cross-section dependence on ^/s. For this paper the reference centre-of-mass 
energy for the weighting is ^/s — 208 GeV instead of = 189 GeV. 

We have also evaluated results based on the number of events passing an optimised cut on 
the Lr value. These results are used as a cross-check, as a direct comparison between data and 
background and as a simpler basis for combining this analysis with other experiments or other 
analyses. 

6 Constraints on New Physics 

We present limits on the pair production of charged scalar leptons, leptonically decaying charged 
Higgs bosons and charginos that decay to produce a charged lepton and invisible particles. 

The 95% CL upper limit on new particle production at ^/s = 208 GeV obtained by com- 
bining the data at different centre-of-mass energies is calculated at each kinematically allowed 
point on a 1 GeV by 1 GeV grid of m and Am, using the Lr distributions for signal and 
background, the Lr values of the data events, and the efficiency of the general selection at that 
point as input. 

In addition to the Monte Carlo statistical error on the signal efficiency, we assess a 10% 
relative systematic error on the estimated selection efficiency to take into account deficiencies in 
the Monte Carlo generators and the detector simulation (5%), uncertainties in the interpolation 
procedure (5%), effect of tau polarisation in stau decay [28] (5%), fluctuations in the shape of 
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the signal Lr distribution (2%), uncertainties from lepton identification efficiency (2%), detector 
occupancy (1%), trigger efficiency (< 1%), and luminosity measurement (0.5%). 

At high values of Am the dominant background to the searches for new physics results 
from W"'"W~ production. High statistics Monte Carlo samples for this process are available 
and these samples describe the OPAL data well [29]. In addition to the Monte Carlo statistical 
error, we assess a 5% relative systematic error on the estimated background to take into account 
uncertainties in the shapes of the Lr distributions and reference histograms, in the interpolation 
procedure, and deficiencies in the Monte Carlo detector simulation. At low values of Am the 
dominant background results from two-photon e^e^l^l~ events. The background uncertainty 
at low Am is dominated by the limited Monte Carlo statistics; the uncertainty is typically 
20-80% at low Am. 

We have examined distributions of the cut variables and background enriched samples in 
assessing these systematic errors. The systematic errors on the selection efficiency and the 
estimated background are considered as global errors applicable to each search and for all 
values of m and Am . In setting limits, the Monte Carlo statistical errors and the systematic 
uncertainties on the efficiency and the background expectation are taken into account using 
numerical convolution. 

6.1 Limits on CTg 

Limits on cTg, the production cross-section for new physics processes multiplied by the branching 
ratio squared, are presented in a manner intended to minimise the number of model assump- 
tions. The 95% CL upper limits at ^/s = 208 GeV shown in Figures 5-10 arc obtained 
by combining the data at the various centre-of-mass energies using the assumption that the 
cross-section varies as for sleptons and charged Higgs and (3/s for charginos, where (3 is 
the particle's velocity in units of c. The chosen functional forms are used for simplicity in 
presenting the data and represent an approximation, most importantly for processes in which 
t-channel exchange may be important, namely selectron pair and chargino pair production. In 
these cases the cross-section dependence on centre-of-mass energy is model dependent, depend- 
ing on the mass of the exchanged particles and the couplings of the neutralinos and charginos. 
The selectron Monte Carlo events were generated with /x = —200 GeV and tan/3 = 1.5 using 
SuSYGEN. We have found by varying // and tan/3, over the range 100 < |//| < 1000 GeV and 
1 < tan P < 50, that the above choice gives a conservative estimate of the selection efficiency 
for selectrons. 

Upper limits at 95% CL on the selectron pair cross-section at y/s — 208 GeV multiplied 
by the branching ratio squared for the decay e~^ c^Xi ^t^^'^ shown in Figure 5 as a function 
of selectron mass and lightest ncutralino mass. These limits are applicable to ClCl and e^e^ 
production. The corresponding plots for the smuon and stau pair searches are shown in Figures 6 
and 7, respectively. Note that if the LSP is the effectively massless gravitino, G , then for prompt 
slepton decays to a lepton and a gravitino the experimental signature would be the same as 
that for i^^i^Xi with a massless neutralino. In this case the limits given in Figures 5-7 for 
m^o = may be interpreted as limits on the decay £~ — >■ £~G. 

The upper limit at 95% CL on the chargino pair production cross-section multiplied by the 
branching ratio squared for the decay xf — > i^Ui (2-body decay) is shown in Figure 8. The limit 
has been calculated for the case where the three sneutrino generations are mass degenerate. The 
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upper limit at 95% CL on the chargino pair production cross-section multiplied by the branching 
ratio squared for the decay xf W^x5 ~^ -^^^Xi (3-body decay) is shown in Figure 9. 

The upper limit at 95% CL on the charged Higgs boson pair production cross-section mul- 
tiplied by the branching ratio squared for the decay H"*" r^Ur is shown as a function of mH+ 
as the solid histogram in Figure 10. The branching ratio for the decay H'*' t^Vt may be 
the dominant one for the charged Higgs masses explored with this data-set. The dashed line 
in Figure 10 shows the prediction from HZHA at = 208 GeV for a 100% branching ratio 
for the decay H='= — > t^v-^. With this assumption we set a lower limit at 95% CL on mH+ of 
92.0 GeV. 

6.2 Expected Limits and Consistency with Expectation 

For each search, wc provide a table showing quantitatively the signal efficiencies and the agree- 
ment of the data with the Standard Model background expectations for a number of values of 
m and Am . Tables 3, 4, 5, 6 , 7 and 8 give the values of the following quantities in the searches 
for selectrons, smuons, staus, charginos with two-body decay, charginos with three-body decay 
and charged Higgs, respectively: 

1. The signal efficiency for the general selection at 208 GeV with statistical error (the effi- 
ciencies at lower energies differ by less than 3%). 

2. The 95% CL upper limit on the cross-section multiplied by the branching ratio squared 
at 208 GeV, obtained by combining the data from each centre-of-mass energy. 

3. The expected 95% CL upper limit on the cross-section multiplied by the branching ratio 
squared in the absence of signal, (0-95). This is calculated using an ensemble of 1000 
Monte Carlo experiments to simulate the data. In each Monte Carlo experiment, the 
total number of candidates is drawn from a Poisson distribution with mean equal to the 
number of events expected from the Standard Model. For each candidate, a value of Lr is 
assigned, chosen randomly according to the expected Lr distribution for Standard Model 
processes. The expected limit at a given point in m and Am is the mean value of the 
limit for the ensemble of simulated experiments. 

4. The confidence level for consistency with the Standard Model, calculated as the fraction 
of the simulated experiments for which the upper limit on the cross-section multiplied 
by the branching ratio squared is greater than or equal to the value calculated using the 
actual data. In the absence of signal, a CL of 50% is expected on average^. 

Note that the number of candidate events with high values of relative likehhood varies 
greatly from search to search and with m and Am . For the selectron and smuon searches where 
particular di-lepton identities are required and the kinematics of the lepton from the slepton 
decay arc measured directly (in contrast to events with tau leptons), relatively few events may 
be present at high relative likelihood. For searches where the signal di-lepton identities arc not 

'^Values of 100% correspond to points where there are no candidate events with non-zero Lr in the OPAL 
data. In this case, all toy Monte Carlo experiments will have a value of (T95 equal to or (if there are Monte 
Carlo candidates with non-zero Xr) greater than the value for the data. 
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that dissimilar to W+W i'^i/i v and a range of lepton momenta are expected, hundreds of 
events may potentially have high values of relative likelihood. 

For some points in m and Am in these tables, the confidence level for consistency with the 
Standard Model is relatively small (of order 1%). The probability of getting a low confidence 
level for one or more points in m and Am for one or more of the search channels depends on 
the degree of correlation among the different (m. Am) points and among the different channels. 
The correlation between adjacent points is strong when the momentum distributions for those 
points are similar. The momentum distributions vary slowly with both m and Am when Am 
is high (hence the clustering of low confidence level values in Table 4), but vary considerably 
with Am when Am is low. 

This effect was investigated by generating 1000 Monte Carlo experiments with the Standard 
Model Monte Carlo as described above. For each experiment and each point on the (m. Am) 
grid at which signal Monte Carlo has been generated, the CL was calculated. It was found 
that 56% of the experiments had at least one point with a CL of 0.9% or less; this is the lowest 
value observed in Tables 3 to 8. 

6.3 Limits on Slepton Masses 

We can use our data to set limits on the masses of right-handed sleptons^ based on the expected 
right-handed slepton pair production cross-sections and branching ratios. The cross-sections 
were calculated using Susygen at each centre-of-mass energy and take into account initial 
state radiation. In Figure 11 we show limits on right-handed smuons as a function of smuon 
mass and lightest neutralino mass for several assumed vahies of the branching ratio squared 
for /i^ /^^Xi- The expected limit, calculated using Standard Model Monte Carlo only, is 
also shown for a branching ratio of 100%. For a branching ratio /i^ — > /x^x° of 100% and for a 
smuon-neutralino mass difference exceeding 4 GeV, right-handed smuons are excluded at 95% 
CL for masses below 94.0 GeV. The 95% CL upper hmit on the production of right-handed 
f^f" multiplied by the branching ratio squared for — >■ r^x^* is shown in Figure 12. The 
expected limit for a branching ratio of 100% is also shown. For a branching ratio — > t^Xi 
of 100% and for a stau-neutralino mass difference exceeding 8 GeV, right-handed staus are 
excluded at 95% CL for masses below 89.8 GeV. No mixing between fj^ and tr is assumed. 
However, the cross-section ratio (7=+=:- /cr=+ =- at \fs 208 GeV varies between 0.90 and 1.17 
depending only on the mixing angle. Using this information, the limits shown in Figure 12 can 
be applied to any degree of stau mixing by multiplying the predicted cross-section for r^f^ by 
the value of a-+-- /a~+~- corresponding to the mixing angle considered. The two broken lines 
in Figure 12 show the range of possible positions of the line defining the excluded region for a 
branching ratio ff" t^Xi of 100% for any degree of stau mixing. 

For the case of a massless neutralino (or gravitino) and 100% branching ratio, right-handed 
smuons and staus are excluded at 95% CL for masses below 94.3 GeV and 89.8 GeV, respec- 
tively, and fi is excluded at 95% CL for masses below 88.7 GeV, for any degree of stau mixing. 

An alternative approach is to set limits taking into account the predicted cross-section 
and branching ratio for specific choices of the parameters within the Minimal Supersymmetric 

^Thc right-handed slepton is expected to be ligliter than tlie left-handed slepton. The right-handed one 
tends (not generally valid for selectrons) to have a lower pair production cross-section, and so conventionally 
limits are given for this (usually) conservative case. 
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Standard Model (MSSM)^. For < -100 GeV and for two values of tan/? (1.5 and 35), 
Figures 13, 14 and 15 show 95% CL exclusion regions in the (m^±, m^o) plane for right- 
handed selectrons, smuons and staus, respectively. For fi < —100 GeV and tan/j = 1.5, right- 
handed sleptons are excluded at 95% CL as follows: selectrons with masses below 97.5 GeV 

for m~- — m^o > 11 GeV; smuons with masses below 91.0 GeV for m~- — m^o > 3 GeV; and 

e Ai ' // Ai ' 

staus with masses below 85.2 GeV for m~- — m^o > 6 GeV. 

6.4 Search for Unequal Mass Particle Production 

All the search results described above are for the pair production of new particles with equal 

mass. New physics signals with unequal mass, such as the production of e^eR in which the two 
selectrons have different mass and each selectron decays to ex^, can potentially be observed in 
the acoplanar di-lepton event sample. 

We have examined the data for di-electron, di-muon and arbitrary di-lepton identity events 
consistent with unequal mass particle production. The search hypothesis is e+e" XY, with 
subsequent decays, X — > £='=N and Y ^ i'^N, where X and Y are two massive charged particles, 
and N is an invisible particle such as the LSP. The free parameters in the search are the masses 
of particles X, Y and N. The kinematics of producing two particles with different mass constrain 
the lepton momenta to be within well defined ranges depending on the new particle masses. 
The data have been examined by scanning the three particle masses in steps of 10 GeV for the 
three different di-lepton identity hypotheses and requiring the measured lepton momenta to be 
within the ranges specified by kinematics. 

The consistency of the data with the Standard Model expectations is then examined. The 
lowest probability mass hypothesis occurs for di-leptons of arbitrary fiavour with particle mass 
mx = 130 GeV, my = 60 GeV and mN = 20 GeV, where 242 events are observed with 197.3 
expected from Standard Model sources (Poisson probability of 0.7%). From an ensemble of 
1000 Standard Model Monte Carlo experiments it is found that 12.5% of the experiments had 
at least one mass hypothesis with a Poisson probabihty of 0.7% or lower. In general, the 
agreement with expectation is good. In conclusion, we have explored this potential weakness of 
the standard equal mass pair production search and found no significant evidence for unequal 
mass particle production. 



7 Summary and Conclusions 

A selection of di-lepton events with significant missing transverse momentum is performed 
using a data sample with an integrated luminosity of 680 pb~^ collected at e"'"e~ centre-of-mass 
energies ranging from 183 to 208 GeV. The observed number of events, 1317, the dependence on 
centre-of-mass energy and the event properties are consistent with expectations from Standard 
Model processes, dominantly arising from W+W~ production with both W bosons decaying 
le])tonieally. 

^In particular regions of the MSSM parameter space, the branching ratio for £^ ^^Xi can be essentially 
zero as a result of competing cascade decays and so it is not possible to provide general limits on sleptons within 
the MSSM on the basis of this search alone. The predicted cross-sections and branching ratios within the MSSM 
are obtained using Susygen and are calculated with the gauge unification relation, Mi = | tan^ 6wM2. 
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Discrimination techniques are employed to search for the pair production of charged scalar 
leptons, leptonically decaying charged Higgs bosons and charginos that decay to produce a 
charged lepton and invisible particles. No evidence for new phenomena is apparent. Upper 
limits on the production cross-section multiplied by the branching ratio squared for each new 
physics process are presented in a manner intended to minimise the number of model assump- 
tions. 

Assuming a 100% branching ratio for the decay i% — > i^Xi: we exclude at 95% CL: right- 
handed smuons with masses below 94.0 GeV for m~- — m^o > 4 GeV and right-handed staus 

with masses below 89.8 GeV for m~- — m^o > 8 GeV. Right-handed selectrons are excluded 
at 95% CL for masses below 97.5 GeV for m~- — m^o > 11 GeV within the framework of 

C Al 

the MSSM assuming < —100 GeV and tan/3 = 1.5. Charged Higgs bosons are excluded at 
95% CL for masses below 92.0 GeV, assuming a 100% branching ratio for the decay H='= — > t'^Ut-. 
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Am 


m-_ (GeV) 


(GeV) 


45 


50 


55 


60 


65 


70 


75 


80 


85 


90 


94 


99 


103 


1 General Selection Efficiency ( % ) | 


2 


5.0±0.7 


3.8±0.6 


4.5±0.7 


2.2±0.5 


2.0±0.4 


1.6±0.4 


1.2±0.3 


0.8±0.3 


0.6±0.2 


0.5±0.2 


0.5±0.2 


O.litO.l 


O.OitO.O 


2.5 


5.0±0.7 


3.8±0.6 


4.5d=0.7 


2.2±0.5 


2.0±0.4 


1.6±0.4 


1.2±0.3 


0.8±0.3 


0.6±0.2 


0.5±0.2 


0.5±0.2 


O.liO.l 


O.OitO.O 


5 


40.7±1.6 


42.9±1.6 


41.9±1.6 


44.5±1.6 


45.0±1.6 


43.5±1.6 


43.1±1.6 


39.9±1.5 


41.6±1.6 


44.1±1.6 


43.1±1.6 


42.3±1.6 


42.4±1.6 


10 


60.9±1.5 


64.8±1.5 


65.4±1.5 


66.4±1.5 


67.7±1.5 


67.2±1.5 


65.9±1.5 


67.8±1.5 


66.0±1.5 


67.4±1.5 


67.4±1.5 


66.6±1.5 


67.8±1.5 


20 


74.8±1.4 


76.3±1.3 


77.4±1.3 


79.1±1.3 


81.2±1.2 


82.4±1.2 


82.5±1.2 


81.1±1.2 


81.4d=1.2 


81.3±1.2 


82.1±1.2 


83.3±1.2 


80.2±1.3 


m/2 


75.0±1.4 


78.6±1.3 


79.4±1.3 


83.2±1.2 


85.9±1.1 


84.4±1.1 


85.5±1.1 


87.8±1.0 


88.2d=1.0 


87.8±1.0 


88.5±1.0 


88.5±1.0 


89.4±1.0 


m-20 


74.9±1.4 


81.7±1.2 


82.1±1.2 


83.6±1.2 


84.2±1.2 


84.9±1.1 


88.4±1.0 


86.0±1.1 


88.1±1.0 


89.9±1.0 


90.9±0.9 


90.2±0.9 


90.8±0.9 


m-10 


74.6±1.4 


79.3±1.3 


79.8±1.3 


82.0±1.2 


85.2±1.1 


85.2±1.1 


89.1±1.0 


89.2±1.0 


88.7±1.0 


90.7±0.9 


90.1±0.9 


90.3±0.9 


90.5±0.9 


m 


69.5±1.5 


77.3±1.3 


80.8±1.2 


83.1±1.2 


85.0±1.1 


88.2±1.0 


87.3±1.1 


89.2±1.0 


89.5±1.0 


88.5±1.0 


91.1±0.9 


90.4±0.9 


91.1±0.9 



95% CL upper limit on cross-section times BR^jc ^cx*,') (fb) 



2 


82.1 


86.6 


99.4 


139 


113 


160 


294 


412 


1230 


2530 


5020 






2.5 


34.0 


33.7 


36.3 


39.6 


43.5 


35.4 


43.7 


58.6 


85.9 


136 


256 


975 




5 


13.8 


11.6 


12.1 


12.3 


9.6 


15.0 


10.3 


11.9 


12.7 


21.4 


23.2 


36.9 


270 


10 


14.0 


17.2 


12.9 


9.7 


10.2 


10.3 


10.0 


10.2 


10.6 


11.1 


14.9 


23.9 


174 


20 


26.9 


23.4 


21.8 


19.1 


16.6 


21.0 


21.5 


16.5 


13.2 


13.0 


14.1 


22.6 


162 


m/2 


34.4 


35.0 


37.2 


34.6 


34.2 


27.9 


25.7 


31.0 


31.2 


37.5 


40.3 


75.6 


282 


m-20 


39.6 


43.1 


49.7 


45.7 


47.6 


53.7 


60.8 


47.7 


35.2 


30.2 


34.3 


60.1 


243 


m-10 


48.3 


56.3 


56.1 


56.5 


49.4 


48.1 


61.7 


54.8 


42.2 


33.9 


40.3 


52.9 


228 


m 


56.0 


59.3 


61.3 


52.0 


51.2 


46.7 


53.9 


62.9 


39.7 


36.5 


37.9 


54.8 


229 


expected upper limit on cross-section times BFi'^ic ^cxi) (fb) 


2 


78.1 


83.2 


98.5 


122 


143 


192 


342 


474 


1360 


2730 


5280 






2.5 


34.4 


38.1 


39.8 


42.0 


43.4 


46.8 


56.0 


72.8 


98.9 


150 


275 


1010 




5 


15.6 


15.4 


15.5 


15.4 


14.4 


15.0 


15.6 


15.9 


16.4 


20.0 


27.0 


41.1 


274 


10 


18.0 


15.6 


14.4 


13.4 


12.7 


12.6 


12.1 


12.4 


13.2 


16.0 


19.9 


29.4 


188 


20 


29.9 


26.7 


24.1 


22.1 


21.3 


18.6 


17.6 


17.2 


17.7 


19.3 


21.8 


27.3 


171 


m/2 


34.0 


32.9 


32.2 


32.8 


31.9 


32.9 


33.4 


35.3 


38.9 


42.1 


43.7 


45.9 


191 


m-20 


36.0 


37.1 


37.8 


39.1 


40.2 


42.6 


44.8 


48.4 


48.3 


49.1 


50.8 


54.9 


219 


m-10 


38.0 


38.3 


37.2 


38.5 


40.3 


44.6 


46.4 


49.5 


48.8 


47.7 


51.1 


55.0 


218 


m 


38.9 


37.5 


37.9 


38.9 


40.8 


43.0 


47.5 


49.7 


47.7 


48.8 


49.7 


53.9 


216 



CL for consistency with SM (%) 



2 


33.4 


34.1 


39.3 


24.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 






2.5 


42.6 


57.2 


53.4 


46.4 


42.8 


63.1 


100.0 


100.0 


100.0 


100.0 


100.0 


100.0 




5 


57.7 


73.5 


72.8 


68.0 


85.2 


44.5 


87.2 


73.9 


71.2 


35.4 


100.0 


100.0 


100.0 


10 


72.9 


32.2 


56.0 


81.0 


67.8 


66.8 


63.3 


69.0 


68.6 


88.9 


100.0 


66.3 


100.0 


20 


54.3 


59.4 


54.8 


58.4 


69.9 


30.6 


22.8 


47.6 


76.0 


87.4 


95.3 


64.2 


40.1 


m/2 


42.2 


36.2 


26.4 


37.2 


36.1 


63.4 


74.5 


58.4 


67.4 


55.1 


52.5 


5.3 


7.2 


m-20 


32.3 


27.5 


17.5 


28.1 


24.9 


19.6 


14.9 


45.6 


79.5 


91.4 


85.1 


32.3 


27.1 


m-10 


18.8 


8.4 


7.3 


9.4 


20.9 


34.5 


14.2 


32.4 


62.0 


79.6 


70.4 


45.9 


33.6 


m 


9.9 


5.3 


4.4 


16.3 


21.5 


33.7 


30.5 


18.8 


65.6 


75.5 


73.6 


42.0 


34.2 



Am 


m-_ (GcV) 


(GeV) 


45 


50 


55 


60 


65 


70 


75 


80 


85 


90 


94 


99 


103 


General Selection Efficiency ( % ) 




2 


8.9±0.9 


8.5±0.9 


7.0±0.8 


4.3±0.6 


5.1±0.7 


2.7±0.5 


2.0±0.4 


0.5±0.2 


0.6±0.2 


0.2±0.1 


0.3±0.2 


O.litO.l 


O.OitO.O 


2.5 


8.9±0.9 


8.5±0.9 


7.0±0.8 


4.3±0.6 


5.1±0.7 


2.7±0.5 


2.0±0.4 


0.5±0.2 


0.6±0.2 


0.2±0.1 


0.3±0.2 


O.litO.l 


O.OitO.O 


5 


53.1±1.6 


53.3±1.6 


50.7±1.6 


52.0±1.6 


51.3±1.6 


50.6±1.6 


47.0±1.6 


48.2±1.6 


48.4±1.6 


48.9±1.6 


47.4±1.6 


47.4±1.6 


46.4±1.6 


10 


71.1±1.4 


73.4±1.4 


70.6±1.4 


71.8±1.4 


70.1±1.4 


69.6±1.5 


70.3±1.4 


65.9±1.5 


73.3±1.4 


70.9±1.4 


69.9±1.5 


70.7±1.4 


70.0±1.4 


20 


82.6±1.2 


82.7±1.2 


84.4±1.1 


83.8±1.2 


84.9±1.1 


83.4±1.2 


85.7±1.1 


82.1±1.2 


84.7±1.1 


83.2±1.2 


83.3±1.2 


83.4±1.2 


83.5±1.2 


m/2 


84.0±1.2 


85.8±1.1 


86.5±1.1 


87.8±1.0 


90.3±0.9 


88.2±1.0 


89.2±1.0 


91.5±0.9 


92.1±0.9 


91.3±0.9 


91.3±0.9 


92.7±0.8 


92.5±0.8 


m-20 


83.4±1.2 


87.1±1.1 


89.4±1.0 


89.6±1.0 


90.0±0.9 


91.3±0.9 


92.6±0.8 


93.4±0.8 


92.4±0.8 


93.5±0.8 


93.3±0.8 


92.9±0.8 


93.3±0.8 


m-10 


88.6±1.0 


87.9±1.0 


89.2±1.0 


90.0±0.9 


91.8±0.9 


91.3±0.9 


93.0±0.8 


90.1±0.9 


93.2±0.8 


93.3±0.8 


93.5±0.8 


94.2±0.7 


94.6it0.7 


m 


88.6±1.0 


88.8±1.0 


88.1±1.0 


91.3±0.9 


92.1±0.9 


93.3±0.8 


91.6±0.9 


93.5±0.8 


93.6±0.8 


92.8±0.8 


92.1±0.9 


93.5±0.8 


94.5±0.7 


95% CL upper limit on cross-section times BK^iji — > /ixV) (fb) 


2 


49.6 


43.4 


52.6 


64.5 


88.0 


152 


242 


522 


1380 


3410 


4820 






2.5 


28.2 


23.0 


23.6 


26.7 


30.1 


34.3 


39.6 


52.7 


61.1 


108 


206 


584 


8850 


5 


11.9 


12.1 


10.0 


10.1 


10.4 


10.9 


11.6 


10.0 


11.6 


14.9 


20.2 


33.3 


245 


10 


14.4 


11.1 


8.4 


7.8 


7.1 


7.8 


6.9 


8.3 


11.3 


14.6 


21.2 


23.7 


173 


20 


21.2 


16.4 


13.3 


12.4 


12.2 


11.6 


11.3 


10.1 


14.3 


16.4 


19.6 


21.9 


146 


m/2 


25.3 


20.6 


18.9 


28.2 


25.2 


34.6 


29.7 


34.5 


22.1 


26.1 


30.6 


49.6 


200 


m-20 


30.4 


28.1 


31.9 


36.5 


38.7 


49.7 


41.0 


36.1 


35.2 


35.6 


36.0 


53.9 


261 


m-10 


28.5 


35.5 


36.7 


43.0 


45.3 


50.1 


43.9 


34.6 


35.7 


34.9 


35.2 


50.8 


241 


m 


31.7 


33.8 


37.0 


43.8 


44.0 


46.4 


36.2 


37.9 


32.3 


32.3 


34.0 


48.4 


222 


expected upper limit on cross-section times BR'^ip MX'i') (fb) 


2 


44.3 


46.5 


56.1 


67.2 


93.0 


145 


229 


494 


1600 


3580 


4970 






2.5 


25.5 


25.1 


27.0 


28.8 


32.1 


36.0 


43.1 


52.3 


72.8 


124 


231 


601 


8870 


5 


12.9 


12.8 


12.8 


12.1 


12.6 


12.9 


12.9 


13.5 


15.8 


19.7 


25.1 


38.6 


267 


10 


13.3 


12.7 


11.8 


11.5 


11.2 


10.7 


10.6 


11.0 


11.7 


13.6 


17.3 


26.1 


177 


20 


20.8 


19.7 


18.9 


17.7 


17.2 


16.6 


15.9 


15.8 


16.0 


17.6 


20.2 


26.8 


160 


m/2 


22.6 


22.9 


23.4 


24.6 


24.7 


26.8 


27.6 


29.6 


32.6 


36.4 


38.6 


44.5 


196 


m-20 


23.3 


25.1 


25.9 


26.8 


28.6 


29.4 


32.1 


35.2 


37.3 


39.9 


42.2 


47.9 


209 


m-10 


24.0 


24.5 


24.8 


25.8 


27.8 


29.4 


32.9 


35.2 


35.5 


39.0 


40.2 


46.6 


222 


m 


23.8 


23.2 


24.3 


26.0 


27.3 


29.5 


32.0 


34.9 


35.8 


37.6 


40.3 


46.5 


218 


CL for consistency with SM (%) 


2 


29.5 


51.4 


55.8 


53.4 


52.9 


44.5 


41.7 


43.5 


100.0 


100.0 


100.0 






2.5 


32.3 


59.9 


63.4 


61.7 


51.1 


60.6 


59.8 


47.7 


100.0 


100.0 


100.0 


100.0 


100.0 


5 


52.0 


47.7 


76.7 


71.4 


71.8 


70.8 


61.9 


68.7 


66.5 


100.0 


100.0 


100.0 


100.0 


10 


31.9 


57.9 


80.2 


84.9 


91.3 


77.9 


92.0 


71.2 


49.8 


35.9 


19.7 


100.0 


100.0 


20 


40.9 


62.6 


81.9 


83.7 


82.1 


83.0 


80.2 


87.4 


54.5 


51.1 


44.3 


70.4 


75.4 


m/2 


30.6 


54.6 


68.9 


29.1 


41.7 


18.4 


34.7 


26.9 


83.5 


80.0 


69.4 


29.7 


39.2 


m-20 


17.2 


29.3 


23.3 


13.5 


14.5 


4.4 


18.7 


42.6 


49.7 


56.1 


61.9 


30.1 


17.1 


m-10 


25.8 


9.9 


8.1 


3.4 


4.9 


3.0 


15.5 


45.9 


42.4 


56.5 


59.3 


31.6 


29.8 


m 


16.6 


9.6 


7.6 


3.5 


5.6 


6.9 


29.0 


34.7 


55.6 


59.1 


63.2 


36.5 


39.5 



Am 


m-_ (GcV) 


(GeV) 


45 


50 


55 


60 


65 


70 


75 


80 


85 


90 


94 


99 1 103 


General Selection Efficiency ( % ) 




2 


O.litO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


2.5 


O.litO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


O.OitO.O 


5 


11.0±0.4 


10.5±0.4 


9.4it0.4 


8.5it0.4 


7.6it0.4 


7.6it0.4 


6.7it0.4 


6.4it0.3 


5.5it0.3 


4.5it0.3 
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Di-lepton identity 



Figure 1: Di-lepton identities for the data at y/s =189-208 GeV compared with Standard 
Model expectations. The data are shown as the points with error bars (statistical errors only). 
The Standard Model Monte Carlo prediction dominated by 4-fermion processes with genuine 
prompt missing energy and momentum (i'^ui'V) is shown as the lightly shaded histogram 
and the background component, arising mainly from processes with four charged leptons in 
the final state, is shown as the darkly shaded histogram. The Standard Model Monte Carlo 
histograms are normalized to the integrated luminosity of the data. The last six bins correspond 
to "single-lepton events". 
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Figure 2: Distributions of (a) the lepton momentum divided by the beam energy, (b) —qcos6, 
(c) acohnearity (in radians) and (d) momentum asymmetry for the event sample produced 
by the general selection at y/s =189-208 GeV. The data are shown as the points with error 
bars (statistical errors only). The Standard Model Monte Carlo prediction dominated by 4- 
fermion processes with genuine prompt missing energy and momentum (i^ui'V) is shown as 
the lightly shaded histogram and the background component, arising mainly from processes 
with four charged leptons in the final state, is shown as the darkly shaded histogram. The 
Standard Model Monte Carlo histograms are normalized to the integrated luminosity of the 
data. In (a) and (b) there are two entries per event for events containing two identified leptons. 
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Figure 3: Distributions of (a) the total visible energy divided by the centre-of-mass energy, 
(b) the net transverse momentum of the event (c) di-lepton mass and (d) recoil mass to the 
di-lepton system for the event sample produced by the general selection at y/s =189-208 GeV. 
The data are shown as the points with error bars (statistical errors only). The Standard 
Model Monte Carlo prediction dominated by 4-fermion processes with genuine prompt missing 
energy and momentum (i'^u i'V) is shown as the lightly shaded histogram and the background 
component, arising mainly from processes with four charged leptons in the final state, is shown 
as the darkly shaded histogram. The Standard Model Monte Carlo histograms are normalized 
to the integrated luminosity of the data. 
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Figure 4: Distributions of the relative likelihood, Lr, for Standard Model Monte Carlo (shaded 
histogram) , signal (open histogram) and data (points with error bars, statistical errors only), in 
the analysis for staus with a mass of 80 GeV for a stau-neutralino mass difference of 60 GeV for 
the a/s = 206.5 GeV centre-of-mass energy bin. The Standard Model Monte Carlo histogram 
is normalized to the integrated luminosity of the data. 
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Figure 5: Contours of the 95% CL upper limits on the selectron pair cross-section times 
BR^ie — >• exi) at 208 GeV based on combining the 183 - 208 GeV data-sets assuming a 
f3^/s dependence of the cross- section. The kinematically allowed region lies below the dashed 
line. The unshaded region at very low Am is experimentally inaccessible in this search. 
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Figure 6: Contours of the 95% CL upper limits on the smuon pair cross-section times BR^{jl 
/iXi) at 208 GeV based on combining the 183 - 208 GeV data-sets assuming a P^/s dependence 
of the cross-section. The kinematically allowed region lies below the dashed line. The unshaded 
region at very low Am is experimentally inaccessible in this search. 
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Figure 7: Contours of the 95% CL upper limits on the stau pair cross-section times BR^{f 
Tx\) at 208 GeV based on combining the 183 - 208 GeV data-sets assuming a P^/s dependence 
of the cross-section. The kinematically allowed region lies below the dashed line. The unshaded 
region at very low Am is experimentally inaccessible in this search. 
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Figure 8: Contours of the 95% CL upper limits on the chargino pair cross-section times branch- 
ing ratio squared for xt ~^ (2-body decay) at y/s = 208 GeV. The hmits have been calcu- 
lated for the case where the three sneutrino generations are mass degenerate. Only sneutrino 
masses above 35 GeV have been considered given constraints from the Z° line-shape. The limit 
is obtained by combining the 183 - 208 GeV data-sets assuming a [3 / s dependence of the cross- 
section. The kinematically allowed region lies below the dashed line. The unshaded region at 
very low Am is experimentally inaccessible in this search. 
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Figure 9: Contours of the 95% CL upper limits on the chargino pair cross-section times branch- 
ing ratio squared for xf ^^z^X? (3-body decay) at y/s = 208 GeV. The hmit is obtained by 
combining the 183 - 208 GeV data-sets assuming a P/s dependence of the cross-section. The 
kinematically allowed region lies below the dashed line. The unshaded region at very low Am 
is experimentally inaccessible in this search. 
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Figure 10: The solid histogram shows the 95% CL upper hmit on the charged Higgs pair pro- 
duction cross-section times branching ratio squared for the decay H^*" t^Pt at y/s = 208 GeV. 
The hmit is obtained by combining the 183 - 208 GeV data-sets assuming a dependence 
of the cross-section. For comparison, the dashed curve shows the prediction from HZHA at 
Y^i = 208 GeV assuming a 100% branching ratio for the decay — t^u^-. The expected hmit 
calculated from Standard Model Monte Carlo alone is indicated by the dotted line. 
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Figure 11: 95% CL exclusion region for right-handed smuon pair production obtained by com- 
bining the -y/s — 183 - 208 GeV data-sets. The Umits are calculated for several values of the 
branching ratio squared for p,^ — > /i^x? that are indicated in the figure. Otherwise they have 
no supersymmetry model assumptions. The kinematically allowed region lies below the dashed 
line. The expected limit for BE? — 1.0, calculated from Monte Carlo alone, is indicated by the 
dash-dotted line. 
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Figure 12: 95% CL exclusion region for right-handed stau pair production obtained by com- 
bining the ^/s = 183 - 208 GeV data-sets. The hmits are calculated for several values of the 
branching ratio squared for t^Xv The selection efficiency for f^f is calculated for 

the case that the decay f~ — > t~Xi pi'oduces unpolarised r^. Otherwise the limits have no 
supersymmetry model assumptions. The two broken lines adjacent to the limit for BB? = 1.0 
show the region in which this limit can vary if stau mixing occurs (sec text). The kinematically 
allowed region is shown by the dashed line. The expected limit for BR^ = 1.0, calculated from 
Monte Carlo alone, is represented by the indicated broken line. 
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Figure 13: For two values of tan (3 and fi < —100 GeV, 95% CL exclusion regions for right- 
handed selectron pairs within the MSSM, obtained by combining the ^/s = 183 - 208 GeV data- 
sets. The excluded regions are calculated taking into account the predicted branching ratio for 
e^Xi- The gauge unification relation, Mi = |tan^6't^M2, is assumed in calculating the 
MSSM cross-sections and branching ratios. The kinematically allowed region lies below the 
dashed line. 
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Figure 14: For two values of tan (3 and /i < —100 GeV, 95% CL exclusion regions for right- 
lianded smuon pairs within the MSSM, obtained by combining the y/s = 183 - 208 GeV data- 
sets. The excluded regions are calculated taking into account the predicted branching ratio for 
fin ~^ f^^Xi- The gauge unification relation, Mi = | tan^ 6'vi/^2; is assumed in calculating the 
MSSM branching ratios. The kinematically allowed region lies below the dashed line. 
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Figure 15: For two values of tan (3 and fi < —100 GeV, 95% CL exclusion regions for right- 
handed stau pairs within the MSSM, obtained by combining the y/s = 183 - 208 GeV data- 
sets. The excluded regions are calculated taking into account the predicted branching ratio for 
T^Xi- The gauge unification relation, Mi = ^tari^ 6wM2, is assumed in calculating the 
MSSM branching ratios. The selection efficiency for f+f^ is calculated for the case that the 
decay t^Xi produces unpolarised r^. The kinematically allowed region lies below the 

dashed line. 
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